Simultaneous sub-second hyperpolarization of the nuclear and electron spins of 
phosphorus in silicon by optical pumping of exciton transitions 
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We demonstrate a method which can hyperpolarize both the electron and nuclear spins of 31 P 
donors in Si at low field, where both would be essentially unpolarized in equilibrium. It is based 
on the selective ionization of donors in a specific hyperfine state by optically pumping donor bound 
exciton hyperfine transitions, which can be spectrally resolved in 28 Si. Electron and nuclear polar- 
izations of 90% and 76%, respectively, are obtained in less than a second, providing an initialization 
mechanism for qubits based on these spins, and enabling further ESR and NMR studies on dilute 
31 P in 28 Si. 
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Enriched 28 Si is the material of choice for silicon-based 
quantum computing schemes involving electron or nu- 
clear spins [l|, 12, y, 14 , lS| , since the removal of the 29 Si nu- 
clear spin results in very long coherence times 0, [3, 0, 13] • 
Several methods for achieving quantum logic with spin 
states of the shallow neutral donor (D ) 31 P in 28 Si have 
been proposed [J [3, El and the manipulation of electron 
and nuclear spin coherences have been demonstrated Q , 
but unsolved challenges include the measurement of sin- 
gle spins and the initialization, or polarization, of these 
spins. Fortuitously, the isotopic enrichment of 28 Si has 
another dramatic effect: the linewidths of many opti- 
cal transitions are drastically reduced [3, [13, EL EH, EH , 
including those involving 31 P. These narrow transitions 
have bee n prop osed both for measurement of single 
spins 13, 13l.ll4| and for preferentially populating specific 

polarization in silicon has 
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spin states 

Electron and nuclear s pin 
been studied for decades [15j 

I23I . [24]], but the nuclear polarization obtained to date 
has typically been less than a few percent, and requires 
thousands of seconds to establish. Very recently, a 31 P 
nuclear polarization of 68% has been reported 24j in a 
high magnetic field, using a variation of a mechanism 
first proposed in 1959 [17], and demonstrated in InSb in 
1963 [18J], but the time constant was still a relatively long 
150 s. The method demonstrated here works at low mag- 
netic field, and can simultaneously hyperpolarize both 
the electron and nuclear spins of 31 P in less than a sec- 
ond. 

Preliminary attempts to demonstrate this mechanism 
yielded relatively small electron and nuclear polariza- 
tions [25| . due to the fact that until now all samples 
of 28 Si with sufficiently high enrichment to resolve the 
donor bound exciton (D°X) hyperfine transitions were 
p-type, with residual boron acceptor concentrations typ- 



ically ten times the 31 P concentration. At low temper- 
ature all donors were therefore ionized (D + ), precluding 
the observation of D° — > D°X transitions, unless above- 
gap excitation provided photoneutralization. This need 
for above-gap excitation has a strong negative effect on 
the achievable electron and nuclear polarizations, since it 
acts to equalize the populations in the four D° hyperfine 
states. 

The results presented here were made possible by a 
newly grown crystal of 28 Si purposely doped with 31 P. 
This n-type sample allowed us to study the D° — > D°X 
transitions without requiring any above-gap excitation, 
resulting in dramatically larger polarizations. An elec- 
tron polarization of 90% and nuclear polarization of 76% 
were obtained simultaneously in less than a second, in a 
magnetic field and temperature regime where the equi- 
librium electron polarization was only ~2%, and the nu- 
clear polarization ~3x 10~ 3 %. These large and rapidly 
produced hyperpolarizations should be sufficient for the 
initialization of qubits, given the existence of algorithmic 
cooling techniques, which expend modest qubit resources 
to quickly boost relatively high initial polarizations above 
levels sufficient for fault-tolerant operation [2II ]. 

The n-type 28 Si sample central to obtaining these re- 
sults was produced from the same 99.991% enriched 28 Si 
as previous [12|, |l3j, |25j p-type samples, except that dur- 



ing the final floating zone growth of the dislocation-free 
single crystal, 31 P was introduced using a dilute mixture 
of PH3 in Ar carrier gas. The resulting crystal had a gra- 
dient in 31 P concentration, and the present results were 
obtained from a slice containing 1 x 10 14 cm~ 3 boron and 
7 x 10 14 cm~ 3 of 31 P. The sample was a disc with (001) 
faces, 1.5 mm thick and approximately 1cm in diameter, 
mounted in a completely strain free manner in a reflect- 
ing cavity, and immersed in liquid He. Other aspects of 
the apparatus and methods have been described previ- 
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FIG. 1: (Color online) The neutral donor bound exciton tran- 
sition and its low-field Zeeman splittings are shown, along 
with a schematic of the selective optical polarization mech- 
anism. On the left are the splittings of the neutral donor 
ground state (D°) and the bound exciton ground state (D°X). 
In the centre the twelve allowed optical transitions between 
the four D° hyperfine states and the four D°X states are in- 
dicated, numbered from 1 to 12 in order of increasing energy. 
On the right is a simplified schematic of the optical polariza- 
tion mechanism for the pump laser tuned to line 6. 



ously [12j, [13|, |25j . The optical polarization of the pump 
and probe lasers do not have a strong influence on the 
results. 

It should be noted that these ultrahigh resolution D°X 
spectra cannot at present be obtained via emission, or 
PL spectroscopy, the method typically used for stud- 
ies of D°X in semiconductors. The spectral resolution 
needed to resolve the hyperfine splitting is beyond the 
capability of available spectrometers. Thanks to tunable, 
single-frequency lasers, the needed high resolution can be 
achieved in absorption mode, which we implement using 
photoluminescence excitation (PLE) spectroscopy. 

The Zeeman spectroscopy of 31 P in 28 Si in the low field 
regime is outlined in FIG. [TJ On the left is a sketch of the 
splittings of the D° and D°X ground states. The Zeeman 
splitting of D°X is due primarily to the projection of the 
hole angular momentum (mh), since the two electrons 
form a spin singlet. D° is split at zero magnetic field by 
the hyperfine interaction into an electron-nuclear singlet 
and a triplet, separated by 117.53MHz (486.1 neV) )2l\ . 
Under a small magnetic field, these become two hyperfine 
doublets, each separated by approximately one half of the 
zero-field splitting. The D° hyperfine states consist of two 
parallel spin states and ||-||) which are not coupled 
by the hyperfine interaction, and two anti-parallel states 
It-IJ.) and life) which are somewhat mixed by the hyperfine 
interaction. The dipole allowed D° — > D°X transitions 
are shown in the centre of FIG.Q] numbered from 1 to 12 
in order of increasing energy, and consist of six hyperfine- 
split doublets. 



The spectroscopic results are summarized in FIG. [U 
beginning at the bottom with the unpumped spectrum 
of the n-t ype sample. This is very similar to previous 
results [HI, 13, [25| for p-type samples, except that the 
transitions are slightly less well resolved, with a full width 
at half-maximum (FWHM) of 220 neV (54 MHz) as com- 
pared to 150 neV (37 MHz). Whether this small increase 
in linewidth is due to concentration broadening from the 
31 P, or to an unintentional reduction of the isotopic en- 
richment, is not known at this time. The populations 
in the four D° hyperfine states are seen to be essentially 
equal under these conditions, as expected. It is impor- 
tant to note that in order to obtain the spectrum at the 
bottom of FIG. [21 a small amount of above-gap excitation 
from a 1047 nm Nd-YLF laser was required in addition 
to the PLE probe laser. This is not because of a need 
for photoneutralization, as for p-type samples, but rather 
because in the absence of any excitation other than the 
PLE probe laser, the PLE probe laser itself strongly po- 
larizes the D° hyperfine populations. This pumping ef- 
fect of the probe laser results in a weak and distorted 
PLE spectrum, but these saturation effects can be cir- 
cumvented by repopulating the hyperfine states equally 
with a small amount of above-gap excitation from the 
1047nm laser. 

For the other PLE spectra shown in FIG. [21 a much 
stronger pump laser (~ 2.8 Wcm -2 ) is set to the desired 
energy while the PLE probe laser (~ 6.5 x 10~ 2 Wcm -2 ) 
is scanned across the transitions. For the pumped spec- 
tra, no above-gap excitation is used, since now the pump 
laser holds the populations of the hyperfine states essen- 
tially fixed, and any above-gap excitation would merely 
reduce this desired polarization. Unfortunately, with- 
out any above-gap excitation the PLE lines are slightly 
broadened and develop a low energy tail. This results 
from Stark broadening due to the ~ 1 x 10 14 cm -3 of 
ionized boron and phosphorus which are present in the 
absence of above-gap excitation. This broadening could 
be made negligible in a sample with a lower concentration 
of boron. 

The dramatic effects of selective pumping on the D° 
hyperfine populations are immediately apparent in com- 
paring the middle two spectra of FIG. [21 to the unpumped 
spectrum. The results for pumping line 6 and 8 are shown 
since these produce the largest, and opposite, polariza- 
tions. The resulting electronic polarization is extremely 
high, with transitions from D° states having the same 
electron spin as the pumped state almost vanishing from 
the spectra. A large nuclear polarization in the hyper- 
fine doublets from the opposite electron spin states is also 
apparent. Identical results, except of course for the ener- 
gies of the transitions, were obtained when the magnetic 
field was doubled and quadrupled. The top spectrum of 
FIG. [21 shows the result when the pump laser was posi- 
tioned at the half height point on the high energy side of 
line 6, referred to as 6'. The results under the conditions 
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TABLE I: Observed D° populations and net electron and nu- 
clear polarizations derived from the spectra shown in FIG. [2] 
and similar spectra for pumping lines 3, 4, 5, 7, 9 and 10. 



1,149,850 1,149,855 
Laser Energy (\ie\!) 



FIG. 2: (Color online) PLE spectra of the 31 P bound exci- 
ton in an n-type 28 Si sample, revealing electron and nuclear 
polarizations obtained by selective optical pumping. At bot- 
tom is a spectrum without any resonant pumping, showing 
essentially equal populations in the four D° states. The up- 
per spectra show the PLE signal when a strong pump field is 
tuned to either line 6 or line 8, or the half-height point on the 
high energy side of line 6, which is indicated as 6'. 



of FIG. are summarized in TABLE [H which gives the 
observed populations of the four D° hyperfine states ob- 
tained by curve-fitting to the spectra, and the resulting 
net electron and nuclear polarizations, when pumping at 
the peaks of lines 3 through 10 (pumping lines 1, 2, 11 
and 12 give results very similar to lines 5, 6, 7 and 8, 
respectively). Results essentially identical to these were 
obtained when the pump and probe intensities were both 
reduced by a factor of one hundred, indicating that high 
intensity is not needed to obtain high polarization. 

The basic mechanism responsible for these polariza- 
tions is outlined on the right hand side of FIG. [T] for 
the case of pumping line 6. Only D° in the |T-IJ-) hy- 
perfine state are converted into D°X, which decay with 
high probability to ionized donors (D + ) plus free elec- 
trons (e - ), due to the dominance of Auger recombina- 
tion [28| for D°X in silicon. Subsequent electron cap- 
ture may then populate the opposite donor electron spin 
state. The pure nuclear relaxation rate is assumed to be 
negligible and therefore not shown, and the pump rate 
is assumed to be much higher than either the optically 
enhanced electron relaxation rate W, or the cross relax- 
ation rate R. Population is therefore removed directly 
from IT-U-), and also from ||f|~) and ||fr), since these are 
coupled to ||JJ.) via R and W relaxation, and builds up 
in |J,JJ-) to the extent that the effective pump rate exceeds 
W . Note that this is very different from earlier mecha- 
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nisms of nuclear polarization based on either saturating 
ESR transitions 15[ or on making the W and R pro - 
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cesses equilibrate at different temperatures 
These earlier methods all involve the bidirectional cou- 
pling of pairs of states, without or with a Boltzmann 
factor, respectively, while selective optical pumping can 
move population unidirectionally from the pumped state 
to the state having opposite electron spin. 

The mechanism limiting the polarization is suggested 
by several observations. When the pump and probe pow- 
ers are reduced a factor of 100 the polarizations are essen- 
tially unchanged. However when the pump laser is tuned 
off a peak the polarizations are reduced substantially, 
as shown for 6' in FIG. [U and TABLE [fl The achiev- 
able polarization is therefore likely limited by nonselec- 
tive photoionization of D° in all four hyperfine states by 
the pump laser due to the tail of the D° photoionization 
continuum. Thus even larger hyperpolarizations could 
be achieved if the substantial inhomogencous broadening 
present in this sample could be reduced. 

A much more detailed model, incorporating the asym- 
metric inhomogeneous broadening, the homogeneous 
linewidth, power broadening, and a dependence on the 
hole state of the D°X has been constructed, which can 
reproduce the details of the spectra, including the spec- 
tral hole burning evident at the pump laser energy in the 
top three spectra in FIG. [2] A full account of the model 
and the fits to various spectra is beyond the scope of this 
Letter, but one important result is the fact that, inde- 
pendent of which component of a given hyperfine dou- 
blet the pump laser is tuned to, the polarization results 
are dominated by the coupling to the higher energy, or 
anti-parallel, component. 

To be useful for quantum computing and magnetic res- 
onance, spin polarization must occur on a reasonably 
fast timescale. Due to the slow detector used in our 
PLE apparatus, the system cannot capture the electron 
polarization dynamics, but is still capable of following 
the time dependence of the nuclear polarization. FIG. [3] 
shows the transient populations of the four D° hyperfine 
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FIG. 3: (Color online) The transient behaviour of lines reveal- 
ing the four D° populations are shown, with the pump laser at 
line 6 in all cases. Pump and probe are switched on simulta- 
neously at t = 0. The dashed curves show the transients with 
the populations in the fully polarized state, while the solid 
curves show the transients when the populations have been 
fully equilibrated. The middle curve in c shows the transient 
after the fully polarized sample has recovered in the dark for 
20 min. 



to which this assumption can be relied upon in the pres- 
ence of strong optical excitation is unclear. 

The rapidly evolving nuclear polarization we observe 
is sufficiently high to be useful for quantum computing. 
Our results also suggest that even higher polarizations 
will be achieved as better samples of 28 Si, with reduced 
inhomogeneous broadening resulting from both higher 
isotopic enrichment and reduced boron content, become 
available. Selective optical pumping may therefore pro- 
vide a viable initialization scheme for the 31 P nuclear 
spin as a qubit in 28 Si, and should work equally well for 
other shallow donors. It should also make possible NMR 
studies on 31 P in 28 Si, both normal NMR using the hyper- 
polarization, as well as optically detected NMR, reveal- 
ing important new information regarding this potential 
qubit. 

This work is supported by the Natural Sciences and 
Engineering Research Council of Canada (NSERC). 



states with the pump laser on line 6. The transient when 
pump and probe are turned on simultaneously is shown 
for both a fully polarized and a fully depolarized initial 
state. The transient from a polarized initial state is a 
fast step function, but from the fully depolarized state 
FIG. [3h, b and c all show an initial overshoot followed 
by a decay to the steady-state value, due to the transfer 
of the hyperfine populations to the dominant ||JJ.) state. 
FIG. [3J1 shows that the population in the ||JJ.) state, on 
the other hand, has a buildup following the initial fast 
transient, with a ^100 ms time scale. FIG.[3t also shows 
the transient behaviour for line 4 when the fully satu- 
rated sample is allowed to recover in the dark for 20 min. 
The partial recovery of the unpolarizcd population indi- 
cates that the nuclear spin relaxation time under these 
conditions is 35 min. 

We have commented already on the fundamental dif- 
ference between selective optical pumping and previous 
methods used to obtain nuclear hyperpolarization. The 
hyperfine-resolved D°X transitions are also advantageous 
in being able to measure directly the relative populations 
in all four hyperfine states, from which any of the net po- 
larizations or their correlations can be determined. ESR 
measures only the population difference between the | Iff) 
and | Iff) states, and between the |tJ|) and |J.-(J.) states, 
while NMR (which has not yet been feasible for dilute 
31 P in Si) measures the population differences between 
|TlT) and |t|)-}7 an d be tween and MJJ-) states. In ana- 
lyzing their ESR data. lMcCamev et all 3] have assumed 
that the |tf|~) vs. |J,f|~) and |T-I|) vs. ||-||) population ra- 
tios are characterized by the same temperature, and that 
this temperature is the lattice temperature. The extent 



[1 
[2 
[3 

[4 
[5. 

[e; 

[-: 

[8 

[9 
[10 

[11 
[12 

[13 

[14 
[15 
[16 

[17 

[18 

[19 
[20 

[21 
[22 
[23 
[24 

[25' 

[26 

[27 
[28 



Electronic address: thewalt@sfu.ca 

B. E. Kane, Nature 393, 133 (1998). 

D. P. Divincenzo, Fortschr. Phys. 48, 771 (2000). 

R. Vrijen et al, Phys. Rev. A 62, 012306 (2000). 

J. J. L. Morton et al, Nature 455, 1085 (2008). 

T. D. Ladd et al., Phys. Rev. Lett. 89, 017901 (2002). 

A. M. Tyryshkin et al, Phys. Rev. B 68, 193207 (2003). 

A. Tyryshkin et al, J. Chem. Phys. 124, 234508 (2006). 

T. D. Ladd et al, Phys. Rev. B 71, 014401 (2005). 

D. Karaiskaj et al., Phys. Rev. Lett. 86, 6010 (2001). 

D. Karaiskaj et al., Phys. Rev. Lett. 90, 186402 (2003). 

M. L. W. Thewalt et al., Physica B 401-402, 587 (2007). 

A. Yang et al, Phys. Rev. Lett. 97, 227401 (2006). 

M. L. W. Thewalt et al, J. Appl. Phys. 101, 081724 

(2007). 

K.-M. C. Fu et al, Phys. Rev. B 69, 125306 (2004). 

G. Feher and E. A. Gere, Phys. Rev. 103, 501 (1956). 

A. Abragam, J. Combrisson, and I. Solomon, C. R. Acad. 

Sc. Paris 246, 1035 (1958). 

G. Feher, Phys. Rev. Lett. 3, 135 (1959). 

W. G. Clark and G. Feher, Phys. Rev. Lett. 10, 134 

(1963). 

G. Lampel, Phys. Rev. Lett. 20, 491 (1968). 

N. T. Bagraev and L. S. Vlasenko, Sov. Phys. Solid State 
24, 1974 (1982). 

A. S. Verhulst et al, Phys. Rev. B 71, 235206 (2005). 

H. Hayashi et al., Phys. Status Solidi c 3, 4388 (2006). 
L. S. Vlasenko et al, Phys. Status Solidi c 3, 4376 (2006). 
D. R. McCamey et al., Phys. Rev. Lett. 102, 027601 
(2009). 

A. Yang et al, ICPS-29 (Rio de Janeiro, 2008, to be 

published in the AIP Conf. Proc. Series) (2009). 

L. J. Schulman, T. Mor, and Y. Weinstein, Phys. Rev. 

Lett. 94, 120501 (2005). 

G. Feher, Phys. Rev. 114, 1219 (1959). 

W. Schmid, Phys. Status Solidi b 84, 529 (1977). 



